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ABSTRACT 

We use a new method of analysis to determine parameters of cosmological 
gamma-ray bursts (GRBs), assuming that their distribution follows the star-formation 
history of the universe. Spectral evolution is calculated from an external shock 
model for fireball/blast wave evolution, and used to evaluate the measured peak flux, 
duration, and vF v peak photon energy for a GRB source occuring at a given redshift 
and with given values of total energy, baryon- loading and environmental parameters. 
We then fit model distributions of GRB sources to the observed peak flux, duration 
and vF v peak photon energy distributions. We find that the observed width of the 
E p and duration distributions can not be explained by cosmological redshift and 
time dilation effects. Rather, broad distributions of total blastwave energies and bulk 
Lorentz factors are necessary to explain the observed distributions simultaneously 
within the framework of our unifying GRB model. We discuss implications for source 
parameter distributions and determine a range of burst parameters consistent with the 
data. 

Subject headings: cosmology: theory — gamma-rays: bursts — gamma-rays: theory — 
radiation mechanisms: non-thermal 

1. Introduction 

The extragalactic origin of GRBs has been confirmed with the discovery of X-ray, optical 
and radio afterglows of GRBs as a result of the Italian-Dutch Beppo-SAX mission (e.g., Costa 
et al. 1997; van Paradijs et al. 1997; Frail 1998). At present, there are 8 GRBs with redshifts 
obtained from emission line measurements of host galaxy counterparts. These are GRB 980425 
at redshift z = 0.0084 flKulkarni et al. 1998a|) , GRB 970228 at z = 0.695 ( p.jorgovski et al. 
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1999a), GRB 970508 at z = 0.835 ( [Bloom et al. 1998a|) , GRB 970828 at z = 0.958 ( pjorgovsk: 



[1999D , GRB 980703 at z = 0.966 ( Pjorgovski et al. 1998Q , GRB 980613 at z = 1.096 ( Pjorgovski 



et al. 1999b]) , GRB 990123 at z = 1.60 ( |Kelson et al. 1999| ) and GRB 971214 at z = 3.418 



( Kulkarni et al. 1998b| ). The redshift inferred for GRB 980425 depends on the validity of 



the GRB 980425/SN1998bw association (Palama et al. 1998|). The host galaxy redshift of 



GRB 970508 supports the original redshift report obtained through absorption line measurements 
in its fading optical counterpart ( Metzger et al. 1997] ), and is strengthened by the recent 



Beppo-SAX announcement ( [Piro et al. 1999| ) of a redshifted iron fluorescence feature found in the 
spectrum of its fading X-ray counterpart between ~ 2.5 x 10 4 and 6 x 10 s following GRB 970508. 

It is not possible to construct a reliable GRB redshift distribution from the statistics of 
eight GRBs. Except for GRB 980425, however, whose nature remains controversial, the obtained 
redshifts place most GRBs at the cosmological epoch of active star formation. Several models for 
the origin of GRBs, such as the collapsar/hypernova scenario (Woosley 1993; Paczyhski 1998) 
or the supranova scenario of Vietri &; Stella (1998), suggest that GRBs are physically related to 
recent star formation. Their cosmological distribution should therefore trace the star formation 
history of the universe ( [Totani 1997| ). An origin of GRBs involving stellar collapse events is also 
in accord with the small measured offsets of fading transient GRB counterparts with respect to 
the disks of the candidate host galaxies ( Kulkarni 1998 , Bloom et al. 1999| ) . By contrast, a larger 



offset is expected in the compact object coalescence scenario (e.g., Eichler et al. 1989; Narayan, 
Paczyhski, & Piran 1992) and, moreover, the redshift distribution of GRBs in this scenario is not 
required to follow the star formation history of the universe due to the large time delays between 
the formation and merging of compact object binaries. 

Several authors have calculated the GRB peak flux distribution, testing the assumption 
of a GRB rate proportional to the observed star formation rate as given, e. g., by Lilly et al. 
(1996) and Madau et al. (1996). Krumholz et al. (1998) demonstrated that besides a GRB 
spatial distribution following the star formation rate of the universe, a variety of other luminosity 
and redshift distributions of cosmological GRBs are consistent with the observed peak flux 
distribution. Totani (1999) used an empirical correlation between the observed time-integrated 
spectral shape of GRBs and their peak flux, and finds that the resulting peak flux distribution is 
only consistent with the star formation history of the universe if the star formation rate between 
z = and 1 is much flatter (SFR[z = 1] ~ 4 x SFR[z = 0]) than deduced from UV observations 
(SFR[z = 1] ~ 15 x SFR[z = 0]), the former evolution being in agreement with independent 
estimates on the basis of galaxy evolution models ( [Totani et al. 1997| ). He finds similar results for 
models invoking binary neutron star or neutron star - black-hole mergers. In contrast, Wijers et 
al. (1998), assuming a standard-candle luminosity of GRBs, find reasonable agreement between 
their theoretical peak flux distribution and the one deduced from BATSE/PVO observations. 

This controversy indicates that the peak flux distribution data alone do not sufficiently 
constrain GRB parameters, in particular their luminosity and redshift distributions. Recently, 
Kommers et al. (1998) have included very faint, "non-triggered" BATSE GRBs in a peak 
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flux distribution study and demonstrated that even this improved data set cannot confidently 
distinguish between a GRB rate tracing the star formation rate and one tracing the redshift 
distribution of AGNs. 

Previous theoretical GRB peak flux distribution studies have either used simple, non-evolving 
representations of the intrinsic burst spectra, such as a thermal bremsstrahlung spectrum 
(Fcnimore et al. 1993) or a single power-law flKrumholz et al. 1998 ), or have summed over a sample 



of observed, time-integrated spectra ( Fenimore &: Bloom 1995 , Wijers et al. 1998j Kommers et al. 



1998 ). Mallozzi et al. (1996) demonstrated that the assumed intrinsic spectral shape of GRBs has 



a significant influence on the results of theoretical GRB peak flux distribution studies. 

The first attempt to combine peak flux distribution studies with other statistical properties 
was done by Fenimore & Bloom (1995), who deduced a typical distance scale for GRBs from the 
observed effect of time dilation on the burst duration by comparing the brightest and dimmest 
BATSE bursts. They found that in order to explain a time dilation by a factor of 2, the dimmest 
bursts had to be located at z > 6. However, they pointed out that the implied isotropic burst 
luminosity was inconsistent with the observed peak flux distribution, even in the case of a strong 
cosmological evolution of the comoving burst rate and of the luminosity if both evolutions are 
parametrized as power-laws in (1 + z). However, the significance of the observed peak flux - 
duration correlation is highly controversial. Mitrofanov et al. (1993) did not find any evidence 
for time dilation between dim and bright bursts, while Norris et al. (1994) deduced a maximum 
redshift of z ~ 2.25 for the dimmest bursts from time dilation effects. 

Here we demonstrate the importance of considering additional statistical properties to 
constrain the intrinsic and environmental properties of GRB fireballs. These include the duration 
distribution and the distribution of peak energies E p of the time-averaged vF v spectra of GRBs. 



The duration distribution of BATSE GRBs shows a pronounced bimodality (Kouveliotou et al 



1993 ) between short (£50 <^ 0.5 s) and long (£50 ^ 3 s) bursts. This has been interpreted as 
evidence for two physically distinct source populations by Katz & Canel (1996), who found that 
the {V/V max ) distribution of the populations of short and long bursts are significantly different 
from each other. Both populations show evidence for a cosmological origin ({V/V ma x) < 0.5), 
but the long bursts with (V/V max ) = 0.282 ± 0.014 are located at larger distances and/or exhibit 
stronger cosmological evolution than the short bursts with (V/V max ) = 0.385 ± 0.019. This was 
confirmed by Tavani (1998) who found that the long / hard bursts show significantly stronger 
deviations from a uniform distribution in Euclidean space than other subclasses of GRBs. 

The vF v peak energies E p vary from burst to burst, although most are between ~ 50 
keV and 1 MeV (Mallozzi et al. 1997; Strohmayer et al. 1998). Assuming that the range of 
E p intrinsic to the burst sources does not evolve with redshift, Mallozzi et al. (1995) found 
that the E p vs. peak flux distribution of BATSE bursts indicates a maximum redshift range of 
(1 + zi)/(l + zioo) = 1.86^0 24 between bursts with peak fluxes of 1 and 100 photons cm~ 2 s -1 , 
respectively. 
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Previous work on modeling the peak flux distribution of cosmological GRBs has mainly been 
carried out in a model-independent way, using phenomenological representations of the burst 
spectral and temporal properties, without specifying physical parameters of the burst sources. In 
particular, detailed predictions of theoretical blast-wave models for the peak power and temporal 
evolution of GRB spectra have never been used to model the observed statistics of GRBs. In this 
paper, we use a parametric discription of the spectrum and spectral evolution predicted by the 
external shock model for cosmological GRBs ( permer et al. 1999a ) to calculate self-consistently 



the peak flux of a GRB at a given redshift, its duration t^o, and the observed peak energy E p of 
the vF v spectrum. Assuming that GRBs trace the star formation history of the universe, this 
analysis helps to constrain not only the total energy of GRBs, but additional parameters such as 
the initial Lorentz factor (or baryon loading factor) Tq of the GRB blast wave, the total energy 
Eq deposited into the blast wave, the density of the circumburster material (CBM), and the 
equipartition factor q, which parametrizes the magnetic field strength and efficiency of electron 
acceleration in the blast wave. 

In Section 2, we present the model equations for a GRB from a fireball/blast wave which 
is energized, decelerates, and radiates by its interaction with a smooth CBM. We apply BATSE 
triggering criteria to this model in order to extract measured values of peak flux, duration and 
peak energy. A function describing the star formation history of the universe is given. In Section 3, 
we calculate peak flux, tso, and E p distributions for a GRB source population involving ranges of 
baryon loading factors and total fireball energies. Fits to the observed distributions are presented. 
The results are discussed in Section 4, and we summarize with implications for source models of 
GRBs. 



2. Peak Fluxes, Peak Energies, and Burst Durations 

Following the treatment given by Dermer et al. (1999a), the photon number spectrum from a 
GRB located at redshift z and its temporal evolution is parametrized by the expression 

[e ' h > 4*<P L e*m e t? [e/e p (t)\- + (v/S) [e/e p (t)]« ' 1 ' 

where e = E/(m e c 2 ) is the dimensionless photon energy, dj, is the luminosity distance to the burst, 
v and 5 are the asymptotic low-energy and high-energy vF v slopes of the GRB spectrum, and t 
is the time in the observer's frame. Throughout this study we use v = 4/3, while a distribution 
of high-energy slopes N(6) oc <5 -0 ' 5 for 5 > 0.05 is assumed (obviously, the normalization of 
the spectral representation Q diverges for 6 < 0). This distribution is in reasonable agreement 
with the distribution of GRB spectral indices with photon number indices > 2 determined in 
the statistical analysis of high-energy spectra of BATSE GRBS by Preece et al. (1998a). The 
measured distribution is only an approximation to the asymptotic high-energy distribution of 
slopes because the BATSE sensitivity range does not always sample photons with energies e S> e p . 
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The vF v peak photon energy is given by 



ep (t) = 3.0.10-^Kf n fOT °^ <1 1/9 (2) 
PW 1 + z \ x-^-v/ 2 for 1 < x < rl /9 . V ' 



Here, q = y^H ( r /4) £,e ^ s the combined equipartition factor, containing the magnetic field 
equipartition factor the electron equipartition factor £ e , and the shock compression ratio r, 
and To is the initial bulk Lorentz factor or baryon loading factor of the blast wave. The term r\ 
is the power-law index of the surrounding CBM density structure, given by the smoothly varying 
function 

n(x) = no x' 71 . (3) 

The CBM is assumed to be isotropically distributed about the source of the GRB. For simplicity, 
we let T) = in the calculations presented here. The term x = R/Rd is the radial coordinate R in 
units of the deceleration radius Rd, and is related to the observing time t by 

f for < x < 1 

1/(29+,) te i<,<rJ*. (4) 



(2g + 1) i - 2g 

Here is the deceleration time scale of the blast wave in the observer's frame, given by 

1/3 

(5) 



u 1 + z 



cr 8/3 




(3 - V) E 



Air no m p c 2 



and g parametrizes the radiative regime. For a non-radiative (adiabatic) blast wave, g = 3/2 — r], 
while g = 3 — r] describes a fully radiative blast wave. Throughout this study, we assume that the 
blast wave is uncollimated. Beaming can be approximately implemented in the prompt and early 
afterglow phases of a GRB by replacing the total energy Eq with the quantity 4n(dEo/dQ), where 
OEq/BO, is the energy radiated per steradian into the solid angle element whose normal is directed 
within an angle <^ 1/Tq of the direction to the observer. The peak vL v luminosity P p (t) depends 
on the parameters of the model according to the relation 

p , A _ c(2g-3 + r ? )(4^m p c 2 ) 1 / 3 1/3 2/3 r s/3 / x 2 ^ for < x < 1 
p{) ~ 2g(v-i + 5-i)(l + z) 2 n ° \x 2 -^9 forl<x<rj /9 . [) 

Using Eq. (|l|), we determine the peak flux $ p (z) of a GRB, averaged over the BATSE trigger 
time scale At in the 50 - 300 keV energy range, by scanning through t\ and finding the maximum 
of the expression 



!h+At e 2 \ 
At I dt S de ^ e ^ z )\i 



(7) 



where e± = 50/511 and e 2 = 300/511. A burst producing a peak flux has a probability P tr (<l > ) to 
trigger BATSE. We parametrize the trigger efficiencies given by Fishman et al. (1994) by 

Ptr($) = exp[-($ /$) a ], (8) 
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where $o an d a depend on the trigger time scale: 



3>o,i024 = 0.26 photons cm~ 2 s _1 , 01024 = 5.3 

^0,256 = 0.53 photons cm -2 s _1 , «256 = 5.7 (9) 

3*0,64 = 1-05 photons cm _2 s _1 , = 5.5 

The different effective trigger thresholds are mainly due to effects of variable background. In 
addition, the trigger criterion corresponding to At = 1024 ms involves a selection bias toward 
bursts with durations ^ 1 s. Due to this duration bias, shorter bursts with the same peak flux are 
less likely to trigger a GRB telescope because the total number of detected photons during the 
trigger time scale is smaller. Thus the effective integration time is determined in this case by the 
burst duration rather than the trigger time scale. 

For a given peak flux value $ p and a given set of GRB parameters, an approximate limiting 
redshift z maiX (& p ) can be defined at which the GRB produces a peak flux $ p (z max ) = <I>o 
corresponding to the respective trigger time scale. Fig. [l] illustrates how z max of a standard 
GRB varies with Tq, q, and 5 for the 1024 ms trigger criterion. For fixed values of q, the peak 
flux increases with increasing Tq until E p attains values greater than the photon energies of the 
detector triggering range. At larger values of Tq, the emission recorded by a detector is dominated 
by the synchrotron emissivity spectrum F v oc z/ 1 / 3 produced by a distribution of electrons with a 
low-energy cutoff. Because this portion of the spectrum is so hard, the peak flux, and therefore 
z max , declines at larger values of To, as shown in Fig. [[]. 

The Lorentz factor fo of a fireball which produces the maximum peak flux in a detector 
sensitive to photons with energy is given by fo ~ 75 [(1 + z)Ed/ (m e c 2 qn^ 2 )] l / 4: , using eq. (|2|). 
The nominal BATSE triggering range, as noted above, is 50 keV < E^ < 300 keV. For values of 
^ Fo, Zmax is anticorrelated with q because for larger q the vF v peak is shifted towards higher 
frequencies, causing the (Fo, z max ) curves to level off at lower values of Tq. The maximum redshift 
is also very sensitive to changes in 5 for Tq ^ Fo- This is a consequence of the fact that a rapidly 
increasing fraction of the radiated energy is emitted at E > 300 keV as 5 approaches 0. Obviously, 
Zmax is also strongly dependent on the radiative regime g, determining the fraction of the available 
energy which is converted into radiation. 

A differential peak flux distribution is calculated from the relation 

2max(*i) 

{ P + ) H J (1 + zf V(l + z) (1 + zf - fl A (2z + z*) ' 1 ' 

2max(<ti+l) 



( Totani 1999 ), where we choose (fio,^A) = (0.3,0.7) and Hq = 65 km s 1 Mpc . In this 



cosmology the luminosity distance dz, is 



dL = £(l + z)[ dz ' (11) 

Ho i J(l + n z')(l + z') 2 -n A (2z' + z' 2 ) 
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Note that Eq. @ refers to a mono-parametric distribution of burst sources, for which a unique 
relation between peak flux and redshift holds. In the lowest flux bin (<3?o < &p < $1), the 
number of detected bursts will be determined by the trigger probability. Consequently, we use 
Zmax(&o) = oo. ugrb{z) is the comoving burst rate per comoving unit volume. We assume that 
ncRB(z) is proportional to the star formation rate (SFR), for which we use a simple representation 
of the function shown in Madau et al. (1998): 

r 1(T 3 • 10 z+1 for z < 1.1 
SFiJ^HMoyr^Mpc- 3 ] = I .13 for 1.1< * < 2.8 (12) 

I 4.2 • 10-°- 4 ^ +1 ) for z> 2.8. 
The vF v peak energy of the burst spectrum as a function of redshift is evaluated at the time t = t& 
of peak power output of the GRB using eq. @. Thus E p /(m e c 2 ) = 3.0 x lQ- 8 qn /2 T%/(l + z). 

For the duration distribution, we calculate £50 as the duration between the times when 25% 
and 75% of the total photon fluence in the 50 - 300 keV band (from eq. Q) has been received. In 
Fig. m we plot £50 as a function of To, q, and 5 for a standard GRB located at z = 1. For small 
values of Tq, for which E p is below or within the detector energy range, £50 is proportional to the 
deceleration time tj, oc T . For larger Tq, E p is above the detector energy range, and the burst 
duration is determined by the time it takes for the vF v peak of the evolving burst spectrum to 
sweep through the detector energy range. For the parameters assumed here, £50 is only weakly 
dependent on the initial bulk Lorentz factor in the high-Fo limit. Using the asymptotic forms of 
eq. (Q) and realizing that most photons are produced after E p (t) has swept through the detector 
energy range, one can analytically show that £50 oc r( 2s+1 )/( 9+7 '/ 8 ) ^ r^ 1 ^ *- 2 ^ 3 ' 1 if r > r , 
where the last expression holds when r/ = 0. 

The burst duration £50 is independent of the equipartition parameter q for small values of Tq 
and weakly positively correlated with q for high Tq. There is a very strong correlation between 
the burst duration and the high-energy spectral index. As 5 approaches 0, £50 increases drastically 
since the very hard tail of the spectrum produces a "7-ray afterglow" which decays with time 
oc t~ x , where % = [45(1 + 5) — 2]/[2g + 1] for r] = 0. This approaches a t" 1 decay if 6 — > and 
g ~ 1.5. The duration is also correlated with the radiative regime g with t§Q decreasing by an 
approximately constant factor over the entire Tq range considered here as g increases. For typical 
parameters, this downward shift can reach a factor of ~ 10 as g spans the range from 1.5 to 3. 

The E p and t§o distributions are calculated simultaneously while scanning through redshift 
space according to Eq. (10). We caution that our analysis does not take into account any effects 



due to the instrumental noise or the diffuse radiation backgrounds recorded by the BATSE 
detectors. The actual £50 and igo durations of a GRB are expected to be slightly longer than 
measured because the additional background noise will dominate the emissions from a GRB at 
early and late times, particularly for weak GRBs. The actual E p distribution of those GRBs which 
trigger BATSE is, however, not expected to be much different from the measured E p distribution, 
because the background spectrum is subtracted in the spectral analyses which yield E p (see, e.g., 
Mallozzi et al. 1995). 
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3. Comparison with Observed Distributions 

We first applied the formalism described in the previous section to calculate peak flux, E p , and 
duration distributions from cosmological GRBs with a single set of values for the burst parameters. 
The theoretical distributions were compared with the peak flux and duration distributions from 
the Third BATSE GRB catalog as compiled by Meegan et al. (1996), and to the E p distribution 
of Mallozzi et al. (1997). The values of E p were evaluated by Mallozzi et al. (1997) by fitting 
the GRB spectrum from the 4B catalog ( Meegan et al. 1997j ) with the Band function ( |Band et 



al. 1993| ). We use the trigger time scale At = 1024 ms, corresponding to a peak flux threshold 
of <&tri,g ~ 0.2 photons cm -2 s _1 . This is the trigger criterion used to extract the peak flux 
distribution data shown in Meegan et al. (1996) and in our Figs. ||, [|, and ||. We note that our 
analysis is properly applied to data sets which are produced by uniform triggering criteria, but 
that the E p distribution might be biased with respect to the duration and peak flux distributions 
since the E p data can be obtained only for the brighter GRBs. Moreover, the E p analysis uses 16 
energy channel data, whereas the peak fluxes are based on the four energy channel discriminator 
data (R. Mallozzi, private communication, 1998). A uniform selection criterion for the peak flux, 
duration and E p distributions should be considered in future studies. 

Figs. H and |I| show the derived distributions of GRB sources for mono-parametric models of 
the blast wave and surrounding medium. The spread of measured values is thus a consequence of 
distance and cosmological effects and the GRB rate history, assumed throughout this paper to be 



proportional to the star formation rate SFR (Eq. [ 12 1 ) . The dependence of the observables on 
the radiative index g which parametrizes the radiative regime and thus determines the luminosity 
of the GRB, are shown in the individual figures. While the redshift and thus the peak flux 
distributions are very sensitive to slight changes in g, the E p and £50 distributions are rather 
insensitive to the radiative regime. A comparison of Fig. [3| and || demonstrates the sensitivity of 
the resulting distributions on the total energy Eq = 10 52 £52 erg. The figures also show that the 
spread in E p and £50 due to cosmological redshift alone is far less than the width of the observed 
distributions, indicating that they are probably dominated by a spread in the intrinsic GRB 
properties. 

We note that our modeling result is not unique. There is at least one ambiguity in the 
sense that virtually indistinguishable peak flux, E p and duration distributions are obtained when 
varying the CBM density no and the initial bulk Lorentz factor Tq, keeping the product uqTq 
constant. The reason for this ambiguity lies in the fact that the peak flux, E p , and all depend 
only on the product uoTq (see eqs. |7]], 0, and ||). An average density of no ~ 10 2 — 10 5 cm" 3 
seems to be appropriate if GRBs are correlated with star-forming regions. 

As can be seen, the peak flux, duration, and E p distributions measured with BATSE requires 
fireballs with energies ~ 10 52 - 10 53 ergs, baryon-loading factors corresponding to To ~ 200 - 300, 
and an equipartition factor q ~ 10 -3 . The mean redshifts of such GRB sources lie typically near 
z ~ 1. The ability of the external shock model to account for the characteristic duration of GRBs 
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was first noted by Rees and Meszaros (1992). 

Fig. [| shows the results of evaluating a large set (N ~ 2000) of mono-parametric burst 
distributions and subsequently adding these distributions according to a variety of intrinsic 
parameter distributions. Apart from the high-energy slopes 5 (see previous section), the total 
fireball energy £52 and the bulk Lorentz factor Tq, are assumed to be distributed according to 
truncated power-laws, N{E^2) oc E^, N(To) oc Tq 7 . The power-law indices e and 7, and the 
boundaries I^mim To.max are free parameters. Comparing the lower limits on £52 deduced from 
observations of different GRBs, such as GRB 970508 (e.g., Waxman 1997), GRB 980703 (Bloom 
et al. 1998b), and GRB 971214 (Kulkarni et al. 1998b), not to mention GRB 980425 (Kulkarni et 
al. 1998a), we know that the total available energy in the blast wave may well vary from source 
to source by several orders of magnitude. We assume that the values of £52 vary in the range 
£52,111111 = 10~ 4 and £52,max = 100. We find that the resulting peak flux, £ p , and £50 distributions 
are only weakly dependent on the actual values of these boundaries. Bursts with £52 <^ 10~ 3 are 
only detectable at very small redshifts where the number density of burst sources is assumed to be 
small due to the strong evolution of the star- formation rate (see eq. |l2|]). They thus contribute 
very little to the observed distributions. For power-law indices e ^ 1, the number of bursts with 
E52 100 is strongly reduced compared to the lower-energy bursts, so that the contribution of 
bursts with energies beyond this value is of minor importance as well. 

Dermer et al. (1999a, 1999b) argue that there could be a wide distribution in the baryon 
loading factor of GRB fireballs, of which preferentially those sources with vF v peaks between 
~ 50 keV - several MeV, corresponding to a baryon loading factor of To ~ To, are detected due 
to the triggering criteria of currently operating GRB monitors and the limitations of telescopes 
at X-ray and ^> MeV energies. Detectability of burst sources with values of Tq <C Tq is strongly 
reduced because most of the flux will be emitted at energies below the sensitivity threshold of 
BATSE. The same is true, though to a lesser extent, for fireballs which produce blast waves 
with To ^> To. Thus we expect that the peak flux distribution observed from a population of 
burst sources with a variety of Tq values will not strongly depend on the actual range of baryon 
loading factors. However, due to the strong dependence of the peak energy E p on Tq through 
eq. (§), we expect that the observed E p distribution may well be influenced by the shape of the 
ro distribution. The dependence of the burst duration on Fq is weaker (see Fig. §), but still 

o In 

considerable. In terms of the deceleration time t^, we have oc T , so that a range in Tq will 
also lead to a broadening of the £50 distribution for a fixed value of uq. 

A reasonable match of our theoretical distributions simultaneously with the peak flux, the E p 
and the t§Q distributions was achieved for Tq max = 260, n = 100 cm" 3 , g = 1.7, and q = 10~ 3 , 
e ~ 1.5, and 7 ~ 0. A low value of q 10 -3 is in agreement with the results of Chiang & Dermer 
(1999), who argue that such a small equipartition parameter is required in order to prevent rapid 
synchrotron cooling. A much larger value of q would result in a low-energy synchrotron spectrum 
$(e) oc e~ 3 / 2 , inconsistent with the observed hard low-energy asymptotes of BATSE GRB spectra 
([Band et al. 1993 , Crider et al. 1997 , Preece et al. 1998b] ) . To constrain the intrinsic distributions 
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in E52 and To, we fix no, g, q, and I^max to the values quoted above, and let the power-law 
indices e and 7 vary as free parameters. In our fitting procedure we neglect x 2 contributions from 
durations £50 < 0.5 s and from E p values of < 100 keV because, as mentioned earlier, bursts with 
E p < 100 keV tend to be too dim to allow the determination of spectral parameters such as E p . 
We obtain a best fit simultaneously to all three distributions for e = 1.52 and 7 = 0.25. Fig. |5| 
shows the peak flux, E p , £50, and calculated redshift distributions from this burst population. Our 
fits are rather insensitive to the actual distribution of To values in terms of To min 

and the index 7 

of the power-law distribution in Tq. GRBs with Tq 220 are preferentially detected for the value 
of no = 100 cm -3 chosen here, and the contribution of dirtier fireballs to the peak flux, E p , and 
£50 distribution is almost negligible. 

Fig. ^ shows the la and 2a confidence contours for the power-law indices e and 7 of the 
£52 and Tq distributions, respectively, if the remaining relevant parameters are fixed to g = 1.7, 
q = 10~ 3 , no = 100cm -3 , and To^ax = 260. The figure confirms that the actual shape of the 
distribution of baryon loading factors is only weakly constrained: A large range of 7 values 
(—0.2 <J 7 <J 0.7) yields acceptable fits to all three distributions. In contrast, the power-law in E52 
is constrained to a narrow range of indices, 1.49 <^ e <J 1.54. 

The normalization of our model peak flux distribution requires nQRB^) = 4.43-10 -5 SFR/M@, 
which yields a local GRB rate tigrbC-z = 0) = 443 GRBs yr _1 Gpc -3 . Assuming a local galaxy 
number density of 4.8 • 10~ 3 Mpc~ 3 flWijers et al. 1998 ), this is equivalent to a local GRB rate of 



92 Galactic events per Myr. This rate is a factor of 3700 higher than the result of Wijers et al. 
(1998). However, we note that the population assumed here contains many undetected bursts for 
the existence of which, consequently, there is no evidence. 

The faintest detectable bursts in our model distributions are located beyond a redshift of 
z ^ 4. The maximum of the redshift distribution is located around z ^ 1, in good agreement with 
the redshifts from GRBs measured so far. The small number of these GRBs does not yet allow an 
independent statistical analysis. 

It is clear from Fig. |^ that the bimodality of the burst duration cannot be explained with 
the continuous burst source population assumed in our model calculations. The £50 bimodality 
is not an instrumental effect, but indicates that GRBs consist of at least 2 separate classes, 
characterized by different physical parameters. The peak flux distribution of Meegan et al. (1996) 
was constructed on the basis of the 1024 ms trigger criterion, which misses most of the short 
bursts. Furthermore, all GRBs for which redshifts could be measured as a consequence of their 
precise localization by the BeppoSAX satellite, belong to the subclass of long bursts. Thus, there 
is very limited information on the short bursts so far. Therefore we do not attempt to model this 
subclass of GRBs in this paper. With new data expected from the upcoming HETE II mission, 
which might enable us to localize also several short bursts and determine their broadband spectral 
characteristics and redshifts, a global statistical analysis similar to the one presented here might 
become possible. 
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The external shock model qualitatively explains the trend (Tavani 1998) that GRBs with 
harder spectra have, on average, smaller (V/V maiX ) values. Consider the simplification where 
fireballs are characterized only by Tq and Eq. Blast waves with larger values of Tq decelerate 
and radiate more rapidly and emit the bulk of their radiation at shorter wavelengths. They are 
therefore more luminous (see Fig. |l|), shorter and harder. More energetic GRBs have longer 
durations (eq. Q). Tavani's long/hard class of GRBs therefore corresponds to the class of fireballs 
with the largest values of Tq and Eq, and these are the ones that can be detected from the greatest 
distances and therefore exhibit the strongest non-Euclidean effects. By contrast, the long/soft 
class of GRBs represents less luminous fireballs with smaller mean values of Tq and Eq. Only being 
detectable from relatively nearby sites, these would display smaller cosmological effects on the 
peak flux distribution. Analogous reasoning would also apply to the short class of GRBs, though 
the bimodality of this population indicates that other parameters must also be taken into account. 

We note that the external synchrotron shock model, as parametrized by eq. (|]), specifically 
applies to the so-called "Fast-Rise, Exponential Decay" FRED-like burst light curves (see Fishman 
& Meegan 1995; Dermer et al. 1999b| ). While the overall spectral shape and luminosity (which 
determines the peak flux and E p distributions) may be similar for more complicated, spiky light 
curves — which could arise from the interaction of the blast wave with an inhomogeneous medium 
( Dermer &: Mitman 1999 ) — the observed durations might be altered by this effect. However, the 
bimodality of the duration distribution is not reasonably explained by a variety of light curves, 
and probably indicates a bimodality in parameter space (Katz fc Canel 1996, Tavani 1998). 



4. Summary and Conclusions 

We have used the cosmological blast wave model to fit simultaneously the peak flux, E p , 
and t§Q distributions as observed by BATSE, using an analytical representation of the spectral 
evolution predicted by the external shock model. The GRB source distribution is assumed to trace 
the star formation history of the universe. For our standard cosmology (Hq = 65 km s _1 Mpc -1 , 
S7o = 0.3, = 0.7), the peak flux and E p distributions, and the t^o distribution of the population 
of GRBs with durations t§o ^ 0.5 s can be modeled with a source population characterized by 
N(T ) oc T °' 25 , r < 260, N(E 52 ) oc £ 5 ~ 2 L52 , 10~ 4 < E 52 < 100, g = 1.7, and q = 10~ 3 , if the 
density of the circumburster material is assumed as hq = 100 cm' 3 . Equally good fits can be 
found for different values of the CBM density, if the product hq Tq is held constant. Our results 
are fairly sensitive to the radiative regime of the blast wave and to the value of the equipartition 
parameter. 

The widths of the observed statistical distributions can not be explained solely by cosmological 
redshift and time dilation effects, but require a broad distribution of intrinsic parameters, in 
particular of the total blastwave energy E52 and the initial bulk Lorentz factor Tq. Assuming 
that these parameters are distributed according to single power-laws, the E52 distribution is well 
constrained for a fixed set of parameters which gives a good fit to the distributions, while the 
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modeling results are rather insensitive to the detailed shape of the Tq distribution. 

Our model calculation implies a local GRB rate of 443 GRBs yr -1 Gpc -3 or 92 Galactic events 
per Myr. This is a higher rate than obtained in some recent estimates (e.g., Wijers et al. 1998) 
and indicates a better chance of detecting relatively young (<^ 10 4 yr) GRB remnants in nearby 
galaxies. The photoionization signatures of such young GRB remnants recently investigated in 
detail by Perna et al. (1999) might then serve as a critical test for the assumed correlation between 
GRBs and star-forming regions. 

The bimodality of the £50 distribution can only be explained in the context of the present 
model if a bimodality in GRB source parameter space is assumed. 

The simplified analytic form used to represent the observed spectra from cosmological blast 
waves provides a reasonable representation of all three statistical distributions considered in this 
paper, thus lending strong support in favor of the external shock model of GRBs. The width of 
the observed E p distribution may be regarded as evidence for the existence of burst sources with 
a wide range of baryon loading factors Tq, as suggested by Dermer et al. (1999a, 1999b). Our 
analysis only weakly constrains the existence of a population of dirty fireballs with lower values of 
Tq, because these dirty fireballs are largely undetectable for BATSE due to the rapid decline of the 
observed peak flux in the BATSE photon energy range with decreasing Tq. BATSE is therefore 
insensitive to a dirty fireball population except for the few events that occur at redshifts z C 1. 
The small number of GRBs with E p ^ 100 keV in the data shown in Figs, pi ||, and || might be, at 
least in part, an artifact of the selection bias mentioned at the beginning of the previous section, 
namely that E p can only be determined for bright bursts, while dirty fireballs, producing spectra 
with low E p , are intrinsically dim. 

Also the existence of a population of clean fireball burst sources cannot be constrained tightly 
on the basis of our results. However, this conclusion is strongly dependent on the radiative regime 
of the average GRB blast wave. For the value g = 1.7 found in our study to be an appropriate 
choice to reproduce the peak flux, E p and duration distributions, an additional population of clean 
fireballs with Tq much larger than the upper limit quoted above (T nax riQ ~ 462 cm -3 / 8 ) can 
not be excluded, although it slightly worsens our modeling results. However, if a radiative regime 
g < 1.6 is assumed, the decline of the peak flux with increasing Tq > Tq becomes much slower, 
implying that clean fireballs would be detected as high-i^ bursts, inconsistent with the observed 
E p distribution. 

We therefore conclude that our analysis is in accord with the existence of a class of dirty 
fireballs with E p <C 100 keV which have generally not been detected with BATSE. Also the 
existence of a significantly larger fraction of cleaner fireballs than detected by BATSE cannot be 
excluded from our analysis, although its consistency with observations depends strongly on the 
best-fit radiative regime of the average GRB blast wave. 

Our results indicate that the most powerful GRBs can be detected at redshifts z p> 4. The 
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slope of ~ 3/2 in the peak flux distribution seems to be a coincidence due to the particular 
distribution of burst source energies, E52, in combination with the non-Euclidean space-time 
geometry. This is in accord with the fact that the burst with the highest redshift measured so 
far, GRB 971214, still has a rather high BATSE peak flux and fluence, while the closest GRB for 
which a redshift estimate exists so far, GRB 980425, was a rather dim burst. 

In summary, we use a parametrization of the blast wave model to calculate the peak flux, 
duration, and vF v peak energy of a GRB with a prescribed set of intrinsic and environmental 
parameters. We assume that the evolution of the GRB rate with redshift is proportional to the 
star-formation history of the universe. By taking into account triggering properties of GRB 
detectors in our calculations, we make a detailed comparison of model distributions with BATSE 
results for the peak flux, t$o, and E p distributions. We have shown how these distributions can 
be self-consistently modeled and used to extract fireball model parameters. Our analysis shows 
that the observed statistical distributions of observables of GRBs are in accord with the scenario 
that GRBs are associated with sites of active star formation, although alternative cosmological 
distributions of GRB sources, which have not been considered in this paper, cannot be excluded. 
Future work must consider evolutionary behaviors to test compact object coalescence scenarios. 
We argue that the short, hard GRBs represent a separate population of burst sources. 

We thank the anonymous referee for very helpful comments which led to considerable 
improvements of the manuscript. We also thank J. Chiang for useful comments. This work is 
partially supported by NASA grant NAG 5-4055. CD acknowleges support from the Office of 
Naval Research. 
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Fig. 1. — The maximum redshift out to which a GRB can be detected above the 1024 ms trigger 
threshold <&u m ~ 0.2 photons cm~ 2 s -1 , as a function of baryon loading factor Tq, equipartition 
factor q, and high-energy vF y spectral index 5. Parameters: Eq = 10 53 erg, no = 100 cm -3 , g = 1.6. 
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Fig. 2. — Burst duration t^o for a standard GRB located at z = 1, as a function of Tq, q, and 



5. Other parameters: Eq = 10 53 erg, no 
deceleration time td for comparison. 



100 cm 3 , g = 1.6. The light solid curve shows the 
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Fig. 3. — Comparison of peak flux, E p k, £50 and redshift distributions for a set of burst parameters 
E52 = 50, Tq = 240, no = 100 cm -3 , q = 5 ■ 10 -4 , for different radiative regimes of the blast 
wave, i. e. different values of g: g = 1.6 (solid histograms); g = 1.8 (dashed histograms); g = 2.1 
(dot-dashed histograms); g = 2.8 (long-dashed histograms). 
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Fig. 4. — Comparison of peak flux, E p k, £50 and redshift distributions for a set of burst parameters 
E 52 = 3.8, T = 220, n = 100 cm" 3 , q = 10~ 3 , 5 = 0.2, for different values of g: g = 1.6 (solid 
histograms); g = 1.8 (dashed histograms); g = 2.1 (dot-dashed histograms); g = 2.8 (long-dashed 
histograms). 
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Fig. 5. — Comparison of our model burst population to the 3B catalog peak flux, E p k, and £50 
distributions, and model redshift distribution. Parameters: 10~ 4 < E52 < 100, e = 1.52, To < 260, 
7 = 0.25, n = 100 cm" 3 g = 1.7, q = 10~ 3 . 
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Fig. 6. — la (solid line) and 2a (dashed line) confidence contours of the parameters e and 7 
parametrizing the distributions of total blast wave energy and baryon loading factor, respectively 
{N(E 52 ) oc E^, N(F Q ) oc To 7 ), for a standard set of values: g = 1.7, q = 10~ 3 , r , max = 260, 
n = 100 cm~ 3 . 



